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QUANTUM CONTINUAL MEASUREMENTS AND A 
POSTERIORI COLLAPSE ON CCR 

V. P. BELAVKIN 


Abstract. A quantum theory for the Markovian dynamics of an open system 
under the unsharp observation which is continuous in time, is developed within 
the CCR stochastic approach. A stochastic classical equation for the posterior 
evolution of quantum continuously observed system is derived and the spon¬ 
taneous collapse (stochastically continuous reduction of the wave packet) is 
described. The quantum Langevin evolution equation is solved for the general 
linear case of a quasi-free Hamiltonian in the initial CCR algebra with a fixed 
output observable field, and the posterior Kalman dynamics coresponding to 
an initial Gaussian state is found. It is shown for an example of the posterior 
dynamics of quantum unstable open system that any mixed state under a com¬ 
plete nondemolition measurement collapses exponentially to a pure Gaussian 
one. 


1. Introduction 

The time evolution of quantum system under an observation which is continu¬ 
ous in time cannot be described by any Schrodinger equation due to the stochastic 
irreversible nature of von Neumann reduction of the wave packet at any instant 
of measurement. An adequate model of the quantum unitary evolution giving a 
continuous collapse by a conditioning with respect to the measurements can be ob¬ 
tained in the framework of quantum stochastic (QS) calculus El, firstly introduced 
for output nondemolition processes in 0EJ and recently developed in a quite gen¬ 
eral form in mil Ej. A stochastic wave equation for an observed quantum system 
derived in J5J by using the quantum filtering method (3j, provides an explanation 
of pure quantum relaxation of an atom under a complete observation 0 (Zeno 
paradox) and a Watch Dog effect [§] f° r the reduced wave function of a quantum 
particle under the continuous observation. 

In this paper we develop a regorous quantum stochastic theory of unsharp nonde¬ 
molition measurements of continual families of arbitrary noncommuting observables 
Rt,x given sequentially in the real space-time (f,x) g R 1+d . In the case d = 0 this 
defines the standard unitary dilation of an instrumental process for the quantum 
measurements, which are continuous in time, considered within an operational ap¬ 
proach by Barchielli and Lupieri |I]. We give the direct proof of stochastic evolution 
equation for the posterior states of a general quantum system under a continuous in¬ 
direct measurement of a noncommutative field-process Rt = {Rt.x | x £ R d }. The 
observed process Y(t) is supposed to be nondemolition in the sense |3] - [£!] of the 
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commutativity [Y (r), X (f)] = 0 of the past observables Y* = {Y(r)\r < t} with the 
Heisenberg operators X(t) of the system for every t and self-nondemolition (com¬ 
mutative) [Y(r), Y(f)] = 0 for all r, t. In that case a posteriori state can be found 
0 for any initial prior state by the Takesaki conditional expectations e{A'(t)|Y t } on 
{Y(r) | r < t}' restricted tothe future von Neumann algebras C t = {X(s)|s > t}". 
We shall show that it is possible to represent the open quantum system under ob¬ 
servation within a class of quantum stochastic evolutions in such a way that the 
observed commutative process Y(t) for the sequential unsharp measurements of a 
noncommutative process Rt is described as the sum of noncommutative Heisenberg 
operators R(t) = U*(t)R t U(t) of the subsystem under the measurement and a clas¬ 
sical (commutative) while noise (error) e(t). The unitary evolution U(t) of such 
systems perturbed by a singular interaction with a meter is described in a ’Bose 
reservoir’ by a quantum stochastic Schrodinger equation DEI, driven by a white 
noise (force) f(t). Note that the force f(t) responsible for the perturbation of the 
system due to the measurements, may appear in the quantum Langevin equation as 
well as the classical (commutative) white noise a-0 But the pair (e, /) cannot 
be described within the classical theory of generalized processes any more because 
the error e(t) does not commute with f(t) given the nondemolition condition for 
R(t) and Y(t) = R(t) + e(t). 

It is interesting to note that stochastic equations of the particular diffusive type 
of 113.31) and <133, in their normalized nonlinear version n m mi, have appeared in 
the physical literature also in connection with phenomenological dynamical theo¬ 
ries of quantum reduction and spontaneous collapse dn - HU- The idea is that 
the wave-function reduction associated to a continual measurement is some kind of 
diffusion process and some particular equations of this type are postulated. Our 
approach shows that this diffusion postulate as well as the continual counting reduc¬ 
tion J3] can be derived in the natural general form from the unitary stochastic evo¬ 
lution of a big quantum system by the conditioning with respect to a chosen nonde¬ 
molition process under the continual measurement. The unsharp self-nondemolition 
measurements and the objectification problem are discussed now intensively in the 
physical literature da m within the Davies-Lewis-Ludwigs operational approach, 
but real progress in clarifying the connection between the operational theory of 
continual measurements [Q and the spontaneous reduction theories m-m can 
be done only by using the quantum stochastic and nonlinear filtering methods 0 
- 0,0 which are considered in this paper. 

2. A QUANTUM STOCHASTIC MODEL WITH CONTINUAL UNSHARP 

MEASUREMENTS 

Let us consider the dynamical problem of a sequential observation in continuous 
time t > 0 of a measurable family L t = {L^ x |x £ A} of operators L x = L t x , x = 
(f,x) in a Hilbert space Ti, where A is a Borel space with a u-algebra A. We do 
not suppose that the operators L x are pairwise commutative or even self-adjoint 
or normal. But we at first assume that they are bounded, L x £ C{TC), almost 
everywhere on the space R+ x A with respect to the product A(da") = dfA(dx) of a 
positive measure A(dx) on the Borel space A and the standard Lebesque measure 
df on R + . Here C(H) denotes the space of continuous (bounded) operators in Ti. 

One can consider for example the problem of the (indirect) measurement of 
spin momenta L(. x = L x , described in the Schrodinger picture by the operators in 
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Ti = C 2 of spin projections L x = \ M(dx)/dx + L*, where M( A) = f A i? x dx is an 
operator-valued measure M( A) £ £(C 2 ) of the momentum in a solid angle A C A 
with R x = L x + L* having the eigen-values ±1, and A(dx) = dx is the standard 
solid angle measure on the sphere A = {x £ R 3 ||x| = 1}, normalized to 4n. 

Due to the absence of a joint spectral resolution for the noncommutative family 
{L x |x £ A}, there is no possibility of measuring the corresponding physical quanti¬ 
ties in the usual (direct) sense. Moreover there is no way within orthodox quantum 
mechanics and measurement theory to describe an observation which is continuous 
in time even for a single self-adjoint operator L x = L with a simple spectrum or to 
predict the dynamics of the quantum system under such an observation due to the 
absence of nontrivial mathematical models for noninstantaneous measurements. 

We shall show that these difficulties can be removed within the quantum the¬ 
ory of open systems and indirect measurements, based on the quantum stochastic 
approach mm The basic idea is that the quantum system under an observation 
must be described as the subsystem of a big system, including a Boson field A as 
a model of an observation channel coupled to the system by a singular interaction. 
The measurement information about the physical quantities L x under such cou¬ 
pling can be continuously extracted in a nondemolition way from the continually- 
sequential unsharp observation of the output field B = UoqAU^ given by the direct 
measurements of the compartible complex observables Z =B+A* 0 . 

Let A(dx) be the Bose-field annihilation measure onTi = R + x A, satisfying the 
canonical commutation relations (CCR) 

( 2 . 1 ) [A(A'),A*(A)]=A(AnA'), VA, A' £ AQA) 

in the Fock space T over the Hilbert space L 2 (T i) of square integrable functions 
of x £ Ti. One can realize Cl T as the space L 2 (T) = ®^L 0 L 2 (T n ) of functions /, 
square integrable in the sense that 

r. n 

I |/(xi,...,X„)| 2 n A(d Xi) < oo, 

0<t n <...<t n Coo 

of chains x = • • • i x n)i %i = (t*,Xj), t\ < ... < t n of all finite lengths |%| = 

n = 0, 1,... with respect to the natural measure A(dy) = ria;e x A(dx). We identify 
the chains % £ T as subsets {aq,..., x n } C Ti, t\ < ... < t n and the time ordered 
elements (x \,..., x n ) £ T” of the n-cube T", so that T = U^L 0 is considered as 
the direct union of the sets T n = {(xi,..., x n )\ti < ... < t n }- Then the annihilation 
operator A(A) of the Boson quanta in a measurable region A £ Ti is 

(2-3) [A(A)f](x)= [ f(x Ux)A(dx), 

J A 

where % LI x is defined as the chain (xi,..., Xj, x,Xj+i,..., x n ) of length n + 1 for 
almost all x = (t , x), namely if t {t\, ..., t n }. 

One can easily find that the operator (IO) is adjoint to the creation operator 
A*(A) of the quanta in A £ -4(Ti) 

(2-4) (A*(A)f)ix) = /(x\ x )’ x(A)=xnA, 

*ex(A) 

with respect to the scalar product ilO) and satisfies the CCR O. where x\x = 
(xi,..., Xj_i, Xi+i,..., x n ) is the complement of the elementary chain x £ Ti in 
the chain % £ T ra with Xj = x £ X- I n fh e following we shall regard the operators 


/ OO /» 

|/(x)| 2 A(dx) = ]T / 

n=O'* 
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-A(A), .A* (A) acting as 12 . 31) . 12.411 in the Hilbert space Tt (g> T of square integrable 
vector-functions h : T —> 74 with the invariant domain T> = U^>i 2?(£), where 

v(0 = {hen®T\ j e w HMx)ll 2 A(dx)<oo}. 

Let us consider the quantum stochastic evolution U t , t £ K+ in H (g> T, given 
by the Hudson-Parthasarathy operator equation [13 Cl dC7 + KUdt = (. LdA * - 
L*dA)U for U(t) = U£ having in our (nonstationary) case the form 

(2.5) dU* + K t U t *dt = [ (dA* (t, dx)L x - L* x dA(t 7 dx))Z7 t *, t/ 0 * = /, 

J A 

where Kt = iH t + A f Lf x L t x \(dx), the integral is taken over x £ A, A(f,E) = 
A([0, t) x E), and dA(f,E) = A(t + d£,E) — A(f,E) is the forward differential of 
the process A(t, E) for fixed E £ A. The necessary condition for the unitarity 
U* = U^ 1 of the family Ut, t > 0 satisfying the quantum stochastic differential 
equation (1231) is m the self-adjointness of the operators H t (Hamiltonian) in Ti 
and that the integrals J A x L t x A(dx) exist and equal K t + K%. 

The solution of equation 12.511 can be described Jj explicitly in terms of the 
quantum stochastic multiple integral in Fock scale provided the conditions 

(2.6) f \\H t \\dt < oo, f j ||L t , x || 2 dtA(dx) < oo, Vs £ R+ 

J t<s J t<s J A 

hold which are sufficient for the existence and uniqueness of the unitary solution 
U t of equation 12.511 with H t = Hf. 

Let us define the output observed process Y (t) of unsharp measurements of the 
continual family {L x + L* | x £ Ti} as the time dependent selfadjoint operator¬ 
valued measure Y(t, E), E £ B on some er-semi-ring B C A with A(E) < oo given 
by the quantum stochastic (forward) differential 

(2.7) dY(t, E) = M (t, E)dt + dQ(t, E), F(0,E) = 0, 

where M(t,E) = f E U t (L t , x + i( jX )f7 t *A(dx), Q(t,E) = A(t, E) + A*(t, E). In the 
case of the initial vacuum state |0) £ T of the Bose field and B generating A, 
the generalized processes Y x {t) = dY(t, dx)/dtA(dx) can be regarded as a com¬ 
plete indirect observation of noncommuting operators R x = L x + L* x , x £ Ti 
given by the instantaneous sequential measurements of the commuting operators 
T(dx) = dY(t, dx). Indeed the differentials dQ(t,E) for all measurable E £ A in 
that case are statisticaly equivalent to Wiener increments with zero mean values 
(0|dQ(£, E)|0) = 0 and minimal covariances (0|dQ(f, E')dQ(£, E)|0) = d£A(E D E') 
compatible with the CCR O- They are independent of the operators M(t, E) = 
f E R x (t) A(dx), defined at the infinitesimal volume E = dx by the Heisenberg oper¬ 
ators R x (t) = U t R x Ut as M(t, dx) = R x (t) A(dx). Hence the differences between 
the increments d Y(t, dx) = Y(t + dt, dx) — Y(t, dx) of the form 12.711 and the op¬ 
erators i? x (£)dfA(dx) are just independent Gaussian variables dQ(t, dx), defining 
the minimal random error of the measurement of the noncommutative family R t = 
{i?t, x |x £ A} in the continuous time t £ R + as white noise Q{t) = {Q x (£)|x £ A}. 
One can consider Y(t, E), E £ B as a coarse-graining Y)(t) = Y(t,Ei) of the fam¬ 
ily Y(t, E), E £ A, corresponding to a a-partition B = {Ej £ A \ i £ 1} of a 
measurable subset M = Ej C A. 
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The following theorem shows that the QS equation (l2~5l) up to the Hamiltonian 
H t corresponds to the unique Evans-Hudson diffusion j(t,X ) = UtXU*, satisfying 
the nondemolition principle 


[A(«),B(i,E)], Vt < s G R+, X € £, E e A 

for all Xft) = jft, X) over a von-Neumann initial subalgebra C C CfH) respectively 
to the given output held 

Bft, E) = f f L x (r) A(dx)dr + Aft, E), E £ A. 

Jo J E 

Theorem 1. Let jft) : C —> CfH (8> T), t £ R+ be a quantum diffusion over a 
unital *-algebra C C CfH) having the QS-differential 

d jft, X) = 7 ft, X)dt + A (t, X)dA*(t) + A *(t, X)dA(t), 

where dA*A = J A dA*(dx)<5 X) A*clA = f A <5 x dAl(dx), and 

Sx(t, X*) = 5* x (t, xy, 7 (f, X*) = 7 (t, xy 

are the linear structural maps C —> £(7i (giA), necessary satisfying the conditions 

(i) + 

(«) K(t,X*X) = j(t,X*)6* x (t,X) + 6 x (t,X*)j(t,X), 

(in) j(t,x*x) = A(t,xyA(t,x)-j(t,xy 7 (t,x)- 7 (t,xyj(t,x), 

with A(X)*A(X) = f A S x (XyS x (X)X(dx), S x (t,I) = 0 = S* x (t,I), 7 (*,/) = 0 . 
The family {X(t) =j(t,X) \ X £ C} satisfies the nondemolition condition 

(2.8) [X(s),Yft,E)} =0, Vs<f 

with respect to the output fields Y(t) £ {B(t, E), B*(t, E)} ; 


dB(t, E) = / j(t, L x )A(dx)dt + dA(t, E), E £ A, 

J E 

dB*(t,Y) = [ jft, L*)A(dx)dt + dA* (t, E), t£R+ 

J E 

if and only if A, A* are the inner differentiations: 

S x (t,X) = j(t,(X,L t ^)), 6*(t,X)=j(t,[Ll x ,X]), 

7 (t,X) = /3(t,X) + ± J j(t,L* x [X,L x ] + [L* x) X]L x ) A(dx), 

where (3(t, X*) = (3(t , A')* zs a jft)-differentiation C —► £(7Y Cg> A): 

P(t, X*X) = j(t, X*)(3(t, X) + (3(t, X*)j(t, A), /3(t, I) = 0. 

In the inner case (3(t, X) = jft, i[H t , A]) these conditions together with \Z.(A) uniquelly 
define the quantum Markov spatial flow j(t,X) = UtXUf given by the Hudson- 
Parthasarathy equatuion IF7> I) . Moreover, the output fields B(t), B*(t) and, hence, 
the nondemolition process Y (t) = B(t) + B*(t) are locally unitary equivalent to the 
input fields Aft), A* ft) and to the commutative process Qft) = {Q(t, E) | E £ B}: 
Bft) = UooA^U^, B*(t) = UooA*(t)U * in the sense 

(2.9) Yft, E) = U s Q(t, E )U;, Vs > t. 

In particular, Y(i,E) = UtY t fE)Uf for all E £ B, where Y t (E) = Q([0,f) x E) = 
Qft, E), Q(d;r) = A(da;) + A*(da;), x £ Ti. 
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Proof. The increments dA'(f) = X(t + dt) — X(t) of the linear *-maps j(t) : X i—> 
X(t), j(t,X)* = j(t,X)* uniquelly define the linear *-map 7 (t) : C —> C(Tt ® T) 
and the adjoint maps A (t), A *(t) due to the linear independence of the differentials 
df and dA*(t), dA(f). An application of the QS Ito formula Ei to the conditions 
j{t,I) = I, j(t,X*X) = j(t,X)*j(t,X) gives 7 (t,I) = 0, A (t,I) = 0 = A*(t,/), 
and 

d (X(t)*X(t)) = dX{t)*dX(t) + d X(t)*X(t) + X(t)*dX(t) 

= [A (t, X)*A(t, X) - 7 (t, X)*j(t , X) - j(t, A)* 7 (t, X)]dt 
+(A*(X)j(X) + j(X)*A*(X))dA 
+ (A (X*)j(X) +j(X*)A(X))dA*. 

Comparing this with the QS differential 

d j(t, X*X) = A*(X*X)dA + A(X*X)dA* - 7 (A*A)d t, 

we obtain the conditions (i), (ii), (iii), found in [?] for the Markovian case. 

If Y(t) is a nondemolition process respectively to X(t), then 

[d X(t), y(a)] = [X(t + dt),Y(s)] - [X(t), F(s)] =0, Vi > s 

for t > s and hence 

[ 7 (t,x),y( s )] = [^(t,x),y( s )] = K(i,x),y( s )] =0, vt > a 

due to the commutativity of di, dA*(i, E), dA(f, E) with Y(s), s < t. Applying the 
QS Ito formula to the condition [X(t),Y(t)] = 0 for Y(t) £ {B(t, E), B*(t, E)} we 
obtain 

d[X(i),B(i,E)] = f{[X(t),L x (t)]-5 x (t,X)) A(dx) = 0 

J E 

d[X(i),B*(i,E)] = f ([X(t),L*(t)] + 5*(t, A))A(dx) = 0, 

J E 

i. e. S x (X) = [X,L x ], S* x {X) = [L*,X] due to [dX(t), Y(t)] = 0, L x (t) = 
j(t, L t , x ), L*(t) = j(t, L* x ). This together with (3(X) = j{i[H , X}) gives y(f, A') = 
j{t,lti x )), where 

7 t (A) = i[H t , X] + i Jj,L* x [X,L x ] + [L*,X]L x ) A(dx) 

is the solution of the equation 

7 t(X*X) — X* r y t (X) — 7t (A)*A = J [L* x ,X][X,L x ] A(dx), 

uniquelly defined up to a ^-differentiation (3 t (X) = i[H t ,X], H t = H t *. The unique 
solution j(t, X) = U t XUt of the derived nonstationary Langevin equation under 
the boundness conditions (ITTHi was found in |7], Corollary 4. 

Let us denote by U(s,t), s>t the solution of the quantum stochastic evolution 
equation (G3J) on the interval [t,s] with U{t,t) = I under the integrability condi¬ 
tions (Troll . The operators U{s,r) commute with Y r , r < s, due to commutativity 
of Y r £ A(r),A*(r) and the operators L t ,L*, d A(t), dA*(t), t £ [r, s) generating 
U(s,r). Hence U S Y t U* = U t Y t U^ because U* = U(s,t)*U* for any s > t and 
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because of unitarity of U(s,t). Using the quantum Ito formula m one can easily 
find 

d Y(t, E) = d{U t Y t {E)U*) = dU t Y t {E)U* t + U t dY t {E)U* t + U t Y t (E)dU* 

+dU t dY t (E)U* + dU t Y t {E)dU t * + U t dY t (E)dU* + dU t dY t (E)dU*, 
d B(t) = dA(t) + U t L t U*d\, dB* (t) = dA* (t) + U t L*U? dX, 

dF(E) = dQ(E)+ [ U(L x + L* x )U*\(dx)=dQ(E) + M(E)dt 
J E 

for Y(t, E) = B{t, E) + B*(t, E), E g B due to the only nonzero infinitesimal 
multiplication dA{t, E')dy4*(t, E) = dtA(EflE / ), where M is defined by 112.71) . The 
relation for the process O is a simple consequence of 12.91) and [X, Q(E)] = 0 
for any E g Al(Ti) and X g C(H) ® Ij? : 

[X(s),Y(t, E)] = [U'XU;, U s Y t (E)U:} = U a [X, r t (E)]££ = 0. 

Remark 1. Considering instead of Y(t) = U s Q{t)U* the sequential measure¬ 
ments of the output momentum process Y(t) = U s V(t)U* s > t, defined by 
V(t,E) = i (A(t,E)-A*(t,E)) as 

dY(t,E) = N{t,E)dt + dV{t,E), Eg B, 

where N{t, E) = 4 f E U(L t , x — L^ X )U^ A(dx), one can extract the information 
about the noncommuting self-adjoint operators S x = 4 ( L x — L*). Moreover, by 
doubling A —> A x {—, +} the space A and considering the family {L t , x -, L tjXj+ } 
with = -bt, x /v7F2 instead of {At, x } one can realize the continuous time- 

sequential indirect observation of the pairs of operators 

R Xt + = —( L x + L x ), R Xt - = ( L x — L x ), x g M_j_ x A 

by the measurement of the two commutative output processes Y^{t) = Uoc Q ^ 

Here Q ^ (t) = A T ( t ) + A \ F (t) are given by the independent Boson measures A T 
on _4(R+ x A) as A T {t,E) = A T ([0,i) x E), Eg B, and U t satisfies the equa¬ 
tion 12.511 with two-fold quantum stochastic integral over A x {—, +} instead of A 
which can be written again as 12.511 in terms of A = {A + + iAJ). The com¬ 

plexified observable process Z = -U ( Y + + iY _) defines the unsharp observation 
Z(t, E) = B(t, E) + A*(f,E), Eg B, A* = (A*j_ + iA*_) of the nondemolition 

output field B(t) = UooAi^U^. 

In the case B = A such the continuous measurement gives a complete nondemo¬ 
lition sequential observation |3] of the non-Hermitian operators L x in terms of the 
complexified output process Z{t) = UooW{t) having the stochastic differential 

(2.10) dZ(t,E)=dt [ L x {t)X{dx) + dW{t, E), Eg B, 

J E 

where L x (t) = U t L tjX Uf and W{t) = (Q + (t) + iQ_(t)) = A(t ) + A*(t) is the 

complex Wiener process in Fock space over L 2 (R + x A x { —, +}) with multiplication 
table 


dW*(t,E)dW(t,E') = dfA(EnE') = dW(f,E')dW*(f,E), 
dW(t,E)dW(t,E') = 0, dW*{t,E)dW*{t,E’) = 0. 
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3. A POSTERIORI QUANTUM DYNAMICS UNDER 
THE CONTINUAL MEASUREMENTS. 

Let us consider the quantum diffusion j(t) : C —> C{H ® T) of the system over 
a unital *-algebra C in H, together with the given nondemolition output fields 
d B = Ld\ + dA, d B* = L*dX + dA*. The operators j{t,X) = X(t) under the 
conditions of Theorem [I] satisfy the quantum Langevin equation 

d X(t) ~dtj\ (L*MX(t), L x (t)] + mt), X(t)]L x {t)) A(dx) 

(3.1) = X(t)]dt+ [ (dA*(i,dx)[X(i), L x (t)] + [L* x (t), X(t)]dA(f, dx)) , 

J A 

having the unique solution X(t) = UtXU *, where C/ t *, t £ K.+ are the unitary 
operators defined by the QS equation 12.511 . and 

K(t) = U t K t U*, K*(t) = U t K*U*, L x {t) = U t L x , t U* t , L* x (t) = U t L^ t U* t . 
The equation EB can be obtained from 12.511 by using the QS Ito formula 
d {U t XU?) = d U t XU? + U t XdU* + dU t XdU* 
and the Hudson-Parthasarathy multiplication table m 

dA*(f,E)dA(f,E') = 0, dA(f, E')dA*(t, E) = dfA(E n E'), 

d A{t, E)dA(t, E') = 0, d A* (; t , E)dA* (t, E') = 0, VE, E' £ A 

The a posterior dynamics of the system under the observation m with a given 
initial state </> 0 is the dynamics <p 0 i—> 7r t , t £ K+ of the a posterior state i r* on £, 

giving posterior mean values Xt = nt{X} of X £ C as stochastic functions of the 

trajectories of the observed process Y l = {y(r)|r < t}. According to the a 
posterior state is defined by the conditional expectation e{X}(t) = et{X(t)\Y 1 } on 
the commutant A/j = {y(r)|r < t}' in C(Tt ® X), which contains Y* and X(t) due 
to the nondemolition property 12.811 . By Theoremn]the operators e{X}(t.) £ A/)' 
have in the Schrodinger picture the form 

(3.2) U*e{X}(t)U t = UfetlUtXU^Y^Ut = J ® 7r t {A}, MX £ C 

since f7 t *A/)'f7 t commutes with C{H) ® I. As a map 7r t : C —> Ad t into the Abelian 
algebra Mt = C C(J~) generated by {y r |r < t} on T , the a posterior state 

satisfies a nonlinear stochastic equation , obtained for the first time with respect 
to y ( t) as the quantum filtering equation in 0] G. Here we shall derive a linear 
quantum stochastic equation for a nonnormalized posterior state <p t {X} = p t n t {X}, 
where p t is a positive stochastic functional p t = p(Y *) of Y r = Q(r), r < t. 
Moreover, we shall prove that the stochastic normalization factor p t can be taken 
as the probability density p(v 4 ) of the trajectories v * = {u(r)|r < t} of the observed 
process Y f ' with respect to the standard probability measure v of a Wiener process 
w, represented in the Fock space T as Q with respect to the vacuum state |0) £ T. 
Once the density operator p t = cf> t {I} is found by the solution of the linear posterior 
evolution equation, the density function /&(■«*) = is given by the Segal (duality) 
transformation Q h-> w of the observable process Q l = U^Y t Ut in the Schrodinger 
picture. 

We shall say the nondemolition observation is complete for the quantum diffusion, 
described by the stochastic evolution equation GEL if the subsets E £ B in inn 
generate the u-algebra A. Let us see now in that case the posterior dynamics is 
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not mixing: i r t = T t (p 0 Tf , i. e. it is defined as fg {X} = (ip ™for <p 0 {X] = 
{ip\Xip) by a posterior stochastic propagator T™ \ip£FL^ ip™ = T^ip. We show 
the renormalized propagator F™ = y/pfT™ also satisfies a linear stochastic wave 
equation F + KFdt = LdwF in FC, given in the Fock space representation by the 
operator evolution equation in Ft ® F, 

(3.3) dF t + K t F t dt = [ L t)X F t dF t (dx), F 0 = I, 

J A 

where L t dY t := f A L t , x dY t (dx) = L t dQ(t), (F t (w)ip | F t (w)ip ) = p?. The proof 
is given in Lemma ^ and Lemma 13 in terms of <b t -{bY} = FpXF t . 

Firstly let us note that the wave propagator F t : Ft —> H <g) Xi t as any other 
adapted Wiener functional of Y is defined in the Fock representation F t = F t |0) by 
the generating functional F* = f r F t (x) Tlxex fl’MMck'O coinsiding with the Wick 
symbol (f\F t \f) = F* for g = f + /*, where \f) £ F is the coherent state 

l/Xx) = e" ll/l|2/2 n f( x )> ll/f = f |/(z)| 2 A(dic) 

xex J 

for a / : Fi —> C with ||/|| 2 < 00 , denoted as / 2 = ||/|| 2 , if f* = /. It helps to 
prove the 

Lemma 1. The solution F t of the stochastic equation O satisfies the equiva¬ 
lency condition F t |0) = Up |0) respectively to the vacuum |0) £ T with the unitary 
propagator Up defined by the equation IW5I) . i.e. F t h = Uph for all h = ip ® |0), 
where ip £ FI. t > 0. 

Proof. To this end we note that ,4-measurability coincides with immeasurability 
in this case and the equation for Fp = Up |0) with Fp = I can be simply obtained 
by allowing the right hand side of (\ref{eq:ccrl.5|) to act on the Fock vacuum |0). 
This gives 

(L t dA*{t) - L* t dA(t))Up\0) = (L t dA*(t) + L t dA(t))Up\0) = L t dY t Fp, 

where L t dA*(t) = f A L t , x dA* (t, dx) due to ^-measurability of the map L t : x 1 —> 
L tjX for almost all t, the commutativity of the increments d.B f (E) = dA(t, E) with 
Up and with F t and the annihilation property df? t (E)|0) = 0 = dA(t, E)10) for 
all E £ A. The equation for the £(7Y)-valued symbol F* = (g|F t |0)e 92 / 2 of the 
nonunitary classical stochastic evolution F t defined by (h g \F t h' 0 ) = (ip\F*ip') for all 
h g = ip® e g2 / 2 \g), h' 0 = ip' ® |0), is given by 

(3-4) ^ F‘ + K t Fg = J L t , x F* g g{t, x)A(dx), g = g* £ L 2 (Ti). 

This coinsides with the equation for (g\Up\0)e 92 / 2 = Fpf t having the same form 
as (Id.HU with Fpf t instead of F t and G t = /q g(r)dr, g(t) = f(t) + f*{t) instead 
of Y t = B t + B* t , I?t(E) = Al(i,E), E £ B, and the initial operator F* 0 = I. It 
means that F* t = Fp and Up |0) = F t |0) due to the uniqueness of the solution of 
the equation dioii proved in 0 under the conditions (E3- □ 

Secondly, let us find the QS Langevin equation for the process X g (t) = U t X g Up, 
Xg = igX, where g £ L^fTi) is a £>(R + )(g)£>-measurable square integrable function 
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and e 4 =: e 9 ^ := e 3t is the Wick ordered exponential 


xexp f ( g(r, x)dQ(r, dx) := e“ e a (gt) 

Jo J A 


of the observable yt(g) = /g <?dW = in the Schrodinger picture, corresponding 
to the product e g {\) = Y\ xex g(x), X € T in the Fock representation of e g = e 9 |0). 
Here and below g t £ Lg(ri) denotes the projection of g with gt{x) = 0, if x ^ 
[0, f) x E for every E £ £>, otherwise = g , and a*(g t ) = f gt(x)A*(dx ) = a(g*)*, 
q{g t ) = (a + a*){gt). Taking into account that this exponential is defined by the 
equation de 4 = e 4 < 7 (t)dW t with eo = 1, we can obtain for G(t) = U t XiU* = 
e g (t)X{t) 

dG + (GK + K*G)dt = dt I {L* X GL X + ( L* X G + GL x )g(x)} A(dx) 

J A 

+ / {[L*, G]d_B(dx) + dB* (dx) [G, L x \ 

J A 

+g(x)(dH(dx) +dH*(dx))G} 

using the quantum Ito formula d(eX) = deX + edX + dedX for e g (t) = U t e f fj U^ 
and ED- It helps to write the equation for the vacuum expectation operator 

&g{X} = (0|G(t)|0) = F t elXF t *, $*{/} = P 9t := <&;{/}, Vs > t 

as (0|{dG+ (GK + K*G)dt}(t)\0) = dt(0\{L*GL + (L*G + GL)g}(t)\0), or equiv¬ 
alently 

A • ‘I-:, 1 /', 


&JX} + *t{K;X + XK t } 


(3.5) 


$ 4 {L* YL t , x + (XL tx + L* X)g(t, x)}A(dx). 


The equation 13.511 . with $°{X} = X , defines both the prior quantum Markovian 
dynamics m M 4 : X i—> F t XFj* as M 4 = $g and an operator-valued generating 
functional P 9 = F^egF^ = lim^oo <f> 4 { I} of factorial (normal ordered) moment 
operators 

(0| :Y(x 1 )...Y(x n ): |0) = 5 n P g /6g(x 1 )... Sg(x n ) | 9=0 
for the measurements at t m < t, x m £ A, m = 1,... ,n of generalized derivatives 
Y(x) = Y(dx)/A(dx) = Y x (t ) of the measure Y(di) on Z?(R + ) ®B. It follows from 
the Weyl representation 

(3.6) e*(g) = exp J^{J g(r,x)dY r (dx) ~ \ J 9i r , x) 2 A(dx)dr} = e 9(9t)_9 ‘ /2 

of the Wick exponent e 4 =: e 9 ^ that equation 13.511 defines the characteristic 
operator 0 4 -j>Y} = (0|e i!/(9t) A'(f)|0) of y{g t ) = Jg t (x)Y(dx) = U t q t (g t )UJ: 

0 4 {A} = F t e iqti9t) XFj = e- s ‘ /2 $‘ g {A}. 

Let us denote by v t = {u t (E)|E £ B} a stochastic trajectory i> t (E) : H —► R 
of the process Y t in the Wiener representation vt( E) = Y t ( E,w) and by i > t (g) = 
f 0 gdv, the Wiener integral of g(x), x £ Ti. Now we prove the absolute continuity 
I t {X}(dw) = &t{X}(vt)dy(uj) of the corresponding instrument I t {X}(E), E £ B 
with respect to the standard Wiener measure d/x(w). 
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Lemma 2. The solution of the equation is given by the expectation 


(3.7) 


& g {X} = J^e Vt M-£/ 2 $ t {X}(v t )dp(u) 


of a stochastic map : X i—> <£> t {X}(i>t) as the nonanticipating function <f>“ = 
&t.(vt) of v r , r <t, normalized by a stochastic operator -function P“ = P(ut) and 
the factorial exponent l?OI) of the representation q{g t ) t—> vt(g). 

Proof. Let us take as the spectrum of the commutative field measure Q(dx), 
denoted as u>(dx) in the standard Wiener representation u>(f) = f f(x)w(dx), f £ 
L 2 (ri) and p as the Gaussian probability measure on f l with the correlations 

[ w(A)w(A')dp(cu) = A(A n A') = (0|Q(A)Q(A')|0> 

J n 

induced by the Fock vacuum state. Then u>(gt) = vt(g) as q(gt) = yt(g) for every 


B -measurable g : Ti 


Ftefg one can obtain 


& g {X} 


and due to U t *|0) = F t |0) and the commutativity e f „F t = 


= (0\U t &iXU;\0) = (0|e‘F t *XF t |0) 


4 (v)F? (ui)XF t (w)d p(ui) 


[ e Vt{9) - 9 * /2 F t *(u)XF t (u;)dp{u). 

J n 


Here F t (uj) = Fjf is the solution F t = Ff of the equation 13.311 as the functional of 
Y r (E) = q( l r (E)), r < t, E £ B in the Wiener representation, where l r (E) is the 
indicator of [0,r) x E, and e‘(u;) = is the Wick exponent 13.(ill . Due to 

arbitrariness of B(M + ) ® B measurable g , it defines the posterior map <f> t = $ t (y t ) 
in <E3 as the classical conditional expectation 


(3.8) 


<f>t{X}(u t ) = [ Ff(w)XF t (u))dg(uj\vt) 
Jn 


with respect to the cr-algebra on D, generated by the data iy(E) = u>([0,r) x E), 
r £ [0,f), E £ B. It is given by integrating on Q with the Gaussian conditional 
measure dg(u>\vt) = dp{u>)/dp{vt), where dp{v t ) is the induced Gaussian probabil¬ 
ity measure on the trajectories v* = {v(r)\r < t} = w l \B of the standard Wiener 
measure w(t, E) = tn([0, t) x E) on B 9 E. Hence the probability measure of the 
data v l for the nondemolition observation E3 with a given initial wave function 
if £ hi has the density 


Pt = 


\Ft{u)if\\ 2 dg(uj\vt) = (if\P (v t )if) = p(v t ), 


where P(u t ) = $t{/}(u) = P“. The non-Gaussian measure d v = pdp defines the 
factorial generating functionals p* = (e*(y)) for the process Y f as (ip\<f>l{I}ip) 
and the mean values (X(t)) of the operators X(t) at the initial states if £ hi as 
(i/'|*h| 0) {A'}i/i) by the averaging 

(if\ & g {X}if) = J e Vt ^- 9 ^ 2 7r t {X}(vt)d^v t ) = cfHX} 
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of the product e‘ {y t )^t{X}{vt), where TT t {X}(v t ) = (i/>\$t{X}(vt)ip)/p(v t ), over 
all the observed in the past trajectories v t . □ 

Let us derive the corresponding linear stochastic equation for the non-normalized 
posterior map iTTSl) X i—> <£> t {X} defining the posterior transformation ^ 0 h^oII 
for any initial <j> 0 by tl t {X} = <b t {X}/p t , p t = ^ 0 {P t }. In the case of a complete 
nondemolition observation it can be obtained in the Schrodinger picture from (ESI) 
in the same way as ED from 12.511 by using the Ito’s formula for FfXF t = $ t {X}: 


d(F t *XF t )+F t (XK t +KfX- / L* tx XL t , x \(dx))F t dt 

J A 


Ft(XL t . x + L*X)F t dY t (dx). 


In the general case the stochastic differential equation for ED gives the following 
theorem. 


Theorem 2. The conditional expectation im defining in El l the absolutely 
continuous operational measure < f , t {X}(ut)d/i(u;) with respect to the Wiener pro¬ 
cess vt(u)), represented in Fock space by Y t = {Yt(E)|E £ B}, satisfies the linear 
stochastic equation 


(3.9) 


d $ t {X} + $ t {XK t + KfX- / L* x XL t , x A(dx)}df 

J A 

/ $ t {XL t}X + L* tx X}dY t (dx) 

J A 


corresponding to the equation ED .for the Wick symbol d>‘{X} = (/|$ t |/), g = 
23?/ t . Lfere L tx are B-measurable operator-valued functions of x £ A, L tx = 0, 
z/x ^ E for any E £ B, defined for almost all t as a conditional averaging of L tx 
with respect to B C A by 

/ Lt, x <?(x)A(dx) = / L* iX 5(x)A(dx) 

Ja ./a 

/or any B-measurable square-integrable g : A i—> R and f(t,x) is defined similary 
by the averaging of f(t, x). In particular, L t x = f E . Lt, x A(dx) /or all x £ Ej, 

if B = {Ej £ _4|i £ /} is a a-partition M = °/ MCA and A(Ej) 7^ 0. 

Proof. By the classical Ito’s formula 

d (e‘$ t {X}) = de‘$ t {A-} + g*d$tW + de*d$ t {A-} 

= J gifix^MX + Xl ttX + L* t x X}dY t (dx) 

-e^tfXKt + K* t X - L* x XL t>x }di 


{XL t}X + L* tx X)g{t , x)A(dx) ^ dt 


we obtain from ED equation ED for & g {X} = (0|e‘$i{X}|0), if we take into 
account that (0|dYt(E)|0) = 0, VE £ B and 

[ (XL ttX + L* t x X)g(t , x)A(dx) = f (XL t . x + L* t x X)g{t, x)A(dx) 
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due to R-measurability of g{t,-). Hence equation (Id.911 describes the conditional 
mean value Id.811 in the Fock representation w* i—> Q* with respect to the probability 
measure dfj,(w l ) induced on M t by the vacuum state: 

J d^w*) = (0|e‘4 t {X}|0). □ 

Remark 2. In the case of the output momentum process, described in the 
Schrodinger picture by Yt(E) = F(t,E),E£B, one can obtain in the same way the 
posterior equation for the non-normalized linear stochastic map <E> t in the form 


d$ t {*} + MXK t + k;x 


L*XL t 


J t,x"^t.xA(dx)}dt 


1 

l 


$ t {xZ t , x - 


L* x A}dY t (dx). 


Then by doubling the space A and considering the time-continuous measurement 
of the commutative family Y t ^ = Qzf(t) as in section ??, one can obtain the 
posterior equation, corresponding to the complex observation Z t ( E) = W(t, E), 
E £ B of Lt = {Tt, x |x £A}: 


d$t{A'} + MXK t + K*X - [ L* x AL t . x A(dx)}df 

J A 

(3.10) = f 4> t {XZ t , x }dZ t *(dx) + f $ t {Z* x X}dZ t (dx). 

J A J A 

In the case B = A of complete complex observation this equation has a factoriz- 
able solution $ 4 {X} = F t *XF t , VAT £ C, where F t , satisfies the stochastic equation 
(EH1) in the complexified version 

d F t + K t F t dt = [ L t ' X F t dZt (dx), Z t = ^= (Y t>+ + iY t ,J). 

J A v 2 


4. A CONTINUAL OBSERVATION OF CCR QUASIFREE DIFFUSION 
Let 5 be a symplectic ((-space, i.e. a complex space with involution 

£>*,)* = £*?«! ,VA,eC 

and skew-symmetric bilinear ((-form s : 5 x 5 —> C, such that s(£^,£) is purely 
imaginary for all £ £ 5: 

s(Z,£,*) = -s(Z\0 ,*(£»£)* = *(£,£»)■ 

We denote by SRS the real space of the vectors £ = £ , by 0 a separating space of 
complex-valued linear functionals $ : £ i—>$(£), on enquiped with the weak* 
topology, 90 = {■$ £ 0 | •d + ^ = 0}, where $**(£) = $(£)*, V£ £ and by 
i?(£),£ £ S an operator {(-representation i?(£)* = i?(£**) of the canonical commuta¬ 
tion relations (CCR) 

(4.1) [i?(£),R(£ # )] = | S (£,£“), V££ 5 

in a Hilbert space 7d associated with a Gaussian state 

(4.2) ^{e^} = = 4(£). 
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Here id o G 30 is defined by the expectation do(C) = ^ol-^O)} °f -R and the 
quadratic form £ 2 = (C,C), satisfying the Heisenberg inequality 

£V - (£, O 2 > \ s{£, rjf , v e R3, 

is defined by the symmetric covariance form 

<£, 0 = i + R(d)R(0} - 0(C)0>(0- 

One can realise i?(£) — $o(C) as double real part 2RHo = (Ho + Hq)(C) of the 
creation operator Hq(£) = H 0 (0)* with the vacuum state </> 0 {X} = (ip 0 \Xi/; 0 ) in 
an initial Fock space H = To over the Hilbert space H = 2* , associated with the 
scalar product 

(CIO = {V, C # ) + \ s{r), C # ), VC, V e 3. 

Indeed, the adjoint operators Ho(0), Hg(C) satisfying the CCR 

[A,(C“MS(0] = m 

generate To by the unitary representation 

(4.3) X(C) = e Wo(4) " ?2/2 e ij4 o (e) e lAo(4) ,C G R2 
of the Weyl operators A'(C) = exp{ii?(C)} on i/j 0 : 

*(COo = MZ)e iA ° ( *Vo, 

X(V)X(0 = e is ^x(n + o, 
and (ip 0 \X(£)i/j 0 ) = 0o(0- 

We shall identify the dual space 0 with the completion of 2 in the (nonde¬ 
generate) norm ICI = (0,C) 1,/2 = \/ (RC) 2 + (3C) 2 on 2, such that $(C) = (C,0- 
Denoting j : C l—> jC = C the canonical bounded map 2 —> H, 

IIJCII 2 = (£10 = I£I 2 + \ «(£, C B ) = I£I 2 + 30 

< ICI 2 + |s(3?C,3C)| < ICI 2 + 2I5RCN3CI < 2ICI 2 , 
one can express the scalar product (CIO through the complex metric bounded oper¬ 
ator g = j*j as (C|0 = (C* 1 , Sv)- Here d = g?/ € 0, V 77 G 2 is the complex functional 
00 = (OlO = (0 Sn) defining together with 0(0 = (C # ,gO* = 010 the #- 
functional 2Rd = 1 ? + 1 ?** = 2 R ?7 + s 3 r 7 , where s : R^, —> 0 is the skew-symmetric 
operator (C,s rf) = s(rj ,C), Isr/I < 2 I?/I due to the Heisenberg inequality. 

Let us consider the quantum diffusion of CCR algebra under the continuous 
measurement of the unbounded operators L x = R(£ x ), x G R+ x A, defined by a 
family {C, x \x G R+ x A} of vectors in 2, weakly square integrable: J* £ T (d^,d)dr < 
00 for all t G R+ and d G 0, where 

(4.4) e t {d\d) = f 0(C t , x )OCf,x)A(dx) = [ |(Cf ,$)| 2 A(dx), 

J A J A 

0(C) = O(RC) + i0(3£) = OC 11 )* f° r & h C € 3. Moreover, we shall suppouse that 
the integral iTOli is a weak* continuous function of 1 ? G 0 such that 

[ /(x)ocS, x )A(dx) = (cf/,o 

J A 
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for every square integrable function / : A —> C, where Ct f is an element of S 
denoted as f A Q x f(x) A(dx). We shall suppose also that the Hamiltonian H t of the 
system under the observation is given in the Fock space Tt by the normal ordering 
H t =: ht{R) : of a quadratic form v t (d) + \ u t (d,d) = h t (d). This means 

= v t (d) + ^ u t {d,d), d = d 0 + 2^t{gij), 

where V , T) = exp{—A ( 77 I 77 ) + Aq(t])}iIi 0 , %/j 0 £ H is the normalized vacuum: A^iji^ = 
0, (ip 0 | ip 0 ) = 1 in the initial space H — J-q- 

Let us suppose that {vt\t £ M+} is a locally integrable family of jj-linear forms 
vt(d) = (d,Vt) and {u>t\t £ M+} is a locally integrable family of real symmetric 
bilinear forms on 30 such that 

\v\\ = / Itvldr < 00 , Ital^ = / luvldr <00 Vt < 00 , 

Jo Jo 

where \vt\ = \Jv%, v\ = {vt,vt), lw t i = sup{wt('# / , d) | \d'\ < 1, li?l < 1}. Assuming 
the weak* continuity of the linear functions v t (d) and ui t (d', d) on 0 9 d, Vd' £ 0, 
we identify Vt(d) with d(vt), Vt £ and uit(d',d), with d'{uitd) = (d',uitd), 
where uit is a symmetric and hence bounded operator on the Hilbert space 0. The 
quadratic form of Ht, corresponding to 

i[H t ,R{Q] = v t {s£) + R(w t s£) , 

gives together with i[R(£),L t , x ] = s(£,Ct,x): A L *t,xi R (0\ = s(d,x^) the linear 
quantum Langevin equation EB for X(t) = j(t, R{0) : 

(4.5) d R(t, £) + R(t, in t st;)dt = d P(t, £) + u*(s£)di. 

Here dP(f,£) = if A {s(£, Ct, x )dA(t, dx) - s(£, C t , x )dA*(t, dx)}, n t : 30 -> fiS is 
the linear imaginary operator n t = A ~f t + iui t , where j t = e t — ej is given by 
the weak* continuous function = 2i^s£ t (d',d) of d,d' £ 30, d'(K t d) = 

Kt(d',d) = (d',n t d), 

Kt(d', d) = i$se t (d', d) + icj t (d' , d), Vd, d' £ 30 . 

The following theorem gives the solution of the operator equation (14.511 together 
with an integral of a S(R.+ ) 0 H-measurable locally squareintegrable function g : 
Ti —* R over the differential 

(4.6) d Y(t, E) = R(t, Ct(E) + Cf (E)) + dQ(E). 

Here £ t (E) £ S , Ct(E) = f E (} x A(dx) = Ct(E) ti is defined for any E for which 
A(E) < 00 due to weak* continuity on 0 9 d of the integral f E d(( t )A(dx). Note, 
that the corresponding unitary quantum stochastic evolution (1.5) with unbounded 
operator 

K t = \ J P(CL)^(C t .x)A(dx) + iHt 

exists only if (ip 0 \K t \ifi 0 ) = A f A ||<^ t x || 2 A(dx) = k t ( 0) < 00 for almost all t. The 
Wick symbol k t (d) = (i/i \K t ip ) is defined in this case as 

k t {d) = k t ( 0) + iv t {d) + ^ (e t (d, d) + iui t (d, d)), Vd = d 0 + 23?(gr?). 
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In this theorem we use the notations I • \\, I ■ l| for the norms 


IkI^ = / l/t r ldr, 1^2 = 


Jo J A 


1/2 


where l£l = (C,^) 1 / 2 for a £ G 9foi, and Ik ( I = IzkjI means the norm Ik^I = 
sup{/I i?I} of the real operator in t on the Hilbert space 0. Let us also de¬ 
note g t {r) = g(r), r <t, g t (r) = 0, / t (r,x) = 0 = / t *(r,x), Vr > t, and 


/*(r,x,£) = g{r,x)+is{£ r ,d x ) = /(r,x,^)*, r < t. 


Theorem 3. Let the equations fTfoTPD for the quantum diffusion on the CCR 
algebra be defined by Vt G IRS, u) t '■ 0 —> IRS and (j x G 0 such that 

\v\\ < oo, lz/«l * < 00, 1C + c^2 < OO, Vt G R+ 

where x (, x £ is a weakly B(R + ) ® B -measurable function of x = (t, x), defined 
by Ct9 = Ct9 f or every g G L|( A). T/ien t/ie equation iO has a unique solution, 
defined in the Hilbert space TL®tF = IFo ® L 2 (T) by the quantum stochastic integral 


(4.7) R(t,0+ V(gt) = f s(v r ,^ r )dr + R(^(t)) + a(ft) + a*(f t ) 

Jo 

along the trajectories £ r = Jf\t,fi), r G [0, t) of the backward predual differential 
equation 

(4.8) -t r +iK r s£ r = / 3 (r,x)(C riX + cJ, x )A(dx), 

Ja 

(4-9) £ t = Ces, c(t) = ^ 9) (t,C) = C 0l 


with g = 0 corresponding to y(gt) = 0. Tfte output integral 

y(gt)= [ [ s(t,x)dy(r,dx) 
do J A 

of a $(R+) x B-measurable function g G Lg(Fi) over f/ie differntial \4--b\ l is given 
also by the quantum stochastic integral 14-ty along the trajectory f r = ipi- 9 ^ ( t, 0) o/ 
f/ie equation with £ = 0, corresponding to R(t,f) = 0. 

Proof. First we write the weak solution of the equation EEJ in the standard 
form 

C r = ^ r (t)C + J ^g(s,x)^ r (s)(C s , x + c!, x )dsA(dx), 

where <p r (i)£ = <Pr°\t,£) is the solution of the equation (14.811 with 5 = 0. The 
resolving operator ip r ( t ) exists as the chronologically ordered exponential 


<Pr(t) = 


K tl S . . . K tn sdti .. . d t n 


n =0 1 


/ r<ti<...<t n <t 


due to the estimate \<p r {t)\ = sup{l</5 r (t)£l | l£l < 1} < 


< 5 >' 


l/« tl l... lK t Jdti... d t n < exp{l2oo — i^\\} 


71—0 


I 0<t\...t n <t 
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because Isl < 2. Hence one can obtain the existence of (£ r ,$) for every r £ [0,i) 
and d £ 0 due to the estimate 

rt 




r J A 


l5(s»x)(¥? r (s)23fJC s ,x^)|dsA(dx) 


< (t)^l + l(?l|IC + C B I|( / l¥>7(s)i?l 2 ds) 1/2 

J r 

< ^l£l + \g\\\Q + C h^y/t — rj Idl exp{l2u> — iyl*}. 

Now we integrate the left hand side in E3, taking into account 14.(ill and 14.XU : 

R{t,0+ [ [ fl'(r,x)(i?(r,23?C r (dx))dr + dQ(r,dx)) 

J o J A 


— R(t , £) + 


23? 


g(r, x)d A(r, dx) > - R(r , — in r sC r )dr 


= R( 0, £ 0 ) (23? {/ g(r, x)d A(r, dx)| + d R(r, j r ) + R(r, iK r s£ r )dr) 

= R(<Po j s) (*,£)) + (23?|^(s(r,x) +*s(^ r ,Cj. iX ))dH(r,dx)| + u r (s£ r )dr), 

where di?(r, £ r ) is the quantum stochastic differential dR(r, C)\^ r , satisfying 14.511 
for t = r. This proves Theorem ©. □ 

Note that the solution R(t,Q of the equation 14.511 given by 14.711 for g = 0 
preserves the CCR G3J and satisfies the nondemolition principle 

(4.10) [R(t, 0, Y g (t)} = 0 , £ 2, g £ 

where Y g (t) = f* f A g(r, x)d Y(r, x). 

It can be proved by using the quantum Ito’s formula and 

[dR(t,0,dR(t,{»)] = [dP(t,t),dP(t,t»)l='y t (Bt,Bt»)dt, 

[d R(t, 0, d Y g (t)} = [d P(t£), dQ ff (i)] = is(£, (Ct + Cl)9(t))dt. 

Indeed, if [_R(t,£), R(t,^)] = \ s(£, £*), then from 14.511 it follows 

[d R(t, Q,R(t, C # )] = nj (s£, s^)dt, [R{t, 0, dl?(^ # )] = -/«t(s£, s^ # )dt, 
and d [R(t, £), R(t, ^)] = (, kJ - n t + 7t)(s£, sC # )df = 0; 

d[R(t,Q,Yg(t)} = [d .R(t,£),Y g (t)} + [i?(t,C),dF s (t)] + [dR(t, £), dY g (t)] 

= [R(t, C),R(t, (C t + Cl)g{t)]dt + is(t c t + cS)ff(*))dt 

= 0 

if [R(t, £), 3C,(i)] = 0 and, hence, [dR(t, £), Y g (t)\ = 0. 

Remark 3. Let 5 = 2“ ® 2 + be an orthogonal decomposition of 2 with respect 
to the complex scalar product , rj) , such that 


(£ ® r?) # = r?* 1 ® £**, s(CKg) = o = s(7? # ,C), VC £ S ,77 £2+. 
spondingly as the negative and pc 

c(C",0 :=w(C",0 >0, VCGS-, 


One can take 2^ correspondingly as the negative and positive subspaces of the 
Hermitian form 2s(3?£, $*£) = is(C,C^) : 
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which is uniquely defined in the case of nondegeneracy: s(£, rf) = 0, V 77 G S => £ = 
0, but it is not obligatory. If ( x G 5“ for almost all 2 ; G Ti, and u(£,\ri) = 0 for 
£ G S - , 77 G S + , then the quantum Langevin equation (14.511 can be written in the 
complex linear form 

(4.11) d L(t, £) + L(t , /c t c£)d t = d A{t, £) + rj t (c£)dt, 

where dA(t,£) = i f A s(£, (l x )dA(t, dx), L(t,£) = R(t,£), c£ = is£, ? 7 t (c£) = 
v t (&£), V£ G 5“ and L(i, £) = 0, dA(f,£) = 0, c£ = 0, ? 7 t (c£) = 0, V£ G S + . The 
equation cna for a complex observation can be written in these terms as 

dZ(t,E) = L(C(E))dt + dW(t,E), E g6. 

5. A CCR QUASI-FREE POSTERIOR DYNAMICS AND CONTINUOUS COLLAPSE 

The solution <m obtained for the quasi-free diffusion 14.5ft with the continuous 
observation 14.(ill of CCR gives the possibility to solve easily the equation 13.51) 
at least for the initial Weyl operators $°{X} = e lR ^ = A(£). To this end let 
us represent the product X (£) ® e* of the operator 14.31) on H = X Q and the 
Wick exponent e‘ = e q<K9t ^~ 9t ^ of the integral yt(g) = q(gt) on T = L 2 (T) as the 
exponent of an operator in Heisenberg picture : 

G(t,£) = e b ’ M X(t,Oe Kgt) = e R{t ^ )+v{9t) -^ /2 , 
where y(g t ) = fg f A g(r,x)dY(r,dx). Due to 14.71) the exponent 

R(t,iQ+y(gt) = s(v r ,(pi 9 \t))dr + R{ip { 0 9) (t)) + a(f*) + a*(/t) j (££) 

can be written in normally ordered form with respect to 

a* (/(*£)) = a*(g - a(f*(i£)) = a{g - 7 s(£ # , V? (ff) (C *£))) 


as 


G(t,Q = c g (t) exp{R(ip { 0 9) {t)) +a(/ t *) + = 

= e ^ (i 5 ) e a * (/tW)) X(i ip[ 9 \t,iC))e aU ^ i ^ ) . 

Here c g (t) = exp{/ 0 ‘ s(v r , £ r d7" - g?/2}, g 2 = f* f A g(r, x) 2 drA(dx) and 

lg(0= lnc g(t) + \ I.M 2 (0> l/t| 2 (0= J q ^/*( 5 x,C)/(r,x,C)drA(dx) 

is given by an integral over the trajectories £ r = ( pi 9 \t , £): 


(5-1) /‘(C) = j Q i S ( V r + Cr) + \ ^ r )}dr , 

where £ \gif) = / A C\^g{t, x)A(dx) = Cfg(i) for every /(-measurable function g(t) : 
xgAh g(t, x). Hence the operator-function d>‘ (£) = <I>g{X(£)}, being the vacuum 
expectation (0|G(f,£)|0) is defined in H by 

$ s(C) = exp{/‘(7C) + -R(<^o ff) (M£))}, 


(5.2) 
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since e a lA)|0) = |0) for every /, where / t *(x) = f*(x) on x £ [0,i) x A and 
ft(r, x) = 0, if r > t. One can easily verify that (15.21) satisfies the equation 13.511 , 
written in the CCR quasi-free case for X = X(£) in the differential form 

+ - («t s£,d) + l - e t (s£,s£)}$‘(£) 

(5.3) = 2^(C t g(t),id+ i s£>$*(£). 

This form of the main equation follows from the relations 

m),xm = s(c,ox(o, \ (x(or(o+r(c)x(0) = <c, wa 

defining the derivative ((,d) of X(£) = e lR ^ and the right hand side in 13.511 by 

R{0*X{0 + X[£)R{0 = ?9f(C, id + \ bZ)X{£), 

i 2 

and also from the definition of the quasi-free Hamiltonian evolution in terms of the 
Weyl operators JH: 

[H,X( 0] = {s(v t ,0 - {u> t s£,id)}X(t) . 

Thus we obtain the solution M*(£) = <!>(,(£) of the Lindblad equation for the CCR 
quasi-free case in term of the characteristic operator M 4 {X(£)} of a prior dynamical 
map (f> 0 i—> having the differential form 15.31) with zero right hand side , as 

well as the operator-valued generating functional P g = $“(0) for the factorial 
moments of the observable process Id.fill 

A posterior quasi-free dynamics of the CCR-algebra under the continual observa¬ 
tion m is described by the characteristic a posteriori function $ t (£) = $ t {X(^)} 
with the Wick symbol (/|$t(OI/) = $g(£)> 9 = f + /* °f the Gaussian form 15.21) . 
Hence the operator-valued function <f> t (£), normalised on the probability density 
operator P t = 4>t(0), in the CCR quasi-free case can be represented as the normal 
ordered functional 15.411 of y t = b t + instead of g = g t , where £ r = E, r {y), in¬ 
stead of £ r = defined by the solution £ r (y) = tp^\t,£) of the backward 

stochastic equation 

(5.4) -d_£ r + in r si r &r = [ (C + Cf. x )dF r (dx), = £. 

J A 

In order to find the operator $ t (£) = $ t (£, y t ) in the form of a function <b t (£, v t ) 
of the trajectories v t (g) = u>(g t ) of the observable process yt{g) = q(gt), let us 
solve the equation with X = X(£), having in the quasi-free case the Wick 

symbol (15.31) in TL = Tq. It can be done in terms of 

<th (£) = (i > v \d = d 0 + 23?(g ? y) , 

by solving the linear stochastic differential equation, coresponding to (15.31) 

*d<M£) + {s(v t ,0 - {Kts£,d) + l - et(s£,s£)}<£ t (£)d t 

(5.5) = j 29f<C tlX ,*3+|0^(0dy t (dx), 

as the equation for a posterior characteristic functuon 0 t (£) = <^{3>t(£)} — fat (0 
with a Gaussian ^> 0 (£) = exp{i$(£) — ^ 2 /2} = 0 lJ (£). The stochastic function (j> t (£) 
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defines the operator-valued function 4> t (£) as the normal ordered form (£): of 
the initial operators R — do = A) + Tig in 7 i. = To- 

Theorem 4. Let the initial state <p 0 of the CCR with a linear quantum 

stochastic evolution 14 have the Gaussian characteristic function 14-2 1 1. Then a 
posteriori nonnormalised state = 0o{^*(£)} = under the continuous 

nondemolition observation m also has a Gaussian form 

(5-6) 0 t (O = p t exp{i$ t (£) - i Pt(€,€)} ■ 

Here 

Pt = ex P J a ^{^(2J?C,x)dF r (dx) - i , dr(23?Cr,x) 2 dpA(dx)} 

is the probability density p t = <j> t (0) = p(yt ) of the observation up to time t, dt(f) = 
(£, fit) is the linear stochastic functional dt = dt(yt) of the posterior mean value of 
R(t , £), satisfying the linear filtering equation 

(5.7) M t ^)+MiK t sOdt= f 2K(e,k t cS )X )dF t (dx)+u t (se)dt 

J A 

with d 0 = d 0 , dit(dx) = dF t (dx) — dt f A d t (29?£ t x ) A(dx), k t = p t +|s, and 
Pt{£,0 = (£,Pt£) is the quadratic form of the posterior covariance of R(t,£), sat¬ 
isfying the Riccati equation with Po(£>£) = (£,£)•' 

(5-8) ^ pt(£, 0 + 2 p t (€, in t si) = e t (s£, s£) - J ^ |25R(£, k t C*, x > | 2 A(dx). 

Proof. Let us find from 15.511 a stochastic equation for A t (i£) = ln<^(£), using 
the Ito’s formula dA t {i£) = (j> t d^> t — \ (dAt(i£)) 2 , and 

(dA t (*C)) 2 = (Cm,) 2 = d t f {<23?C t , x , XM)) + (i% c t # , x , sO} 2 A(dx) 

J A 

= +2K t (A , t (i0,s0 - P(s£,s£)}di 

due to (dYt(dx)) 2 = dtA(dx), where A t (£) = d \ t (£), 

ihCM) = / (2^Ct,x)^) 2 A(dx), iq(s£,s£) = j (9£ t x ,s£) 2 A(dx), 

J A J A 

K t (d,sf) = /(2SRC ti x, , !?)(*^Clx ) s0A(dx) = (k t s£,d). 

J A 

It gives a quasilinear stochastic equation of the first order for A t (£): 

dA t + {s(£, u t ) + (*rc t s£, A*) - ^ £t(s£, s£)}dt 
= / A {(2^,x,A t ) + (3Ct # ,x,sO}dTt(dx) - 1 A t )dt, 
where k t = K t — k t and it = e* — e t . This equation has a quadratic form solution 
At(0 = lnp t + dt(£) + - Pt(£,Q, X t (0 = d t + Pt£, A t = p f , 
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where 

din p t = dA t (0) = - \ M^t,^t)dt, 

dd t = dAj(O) = f (2p t 9?£ t x + s3C tx )dYt(dx) - {(p t /x t -isKj)d t -sv t }dt, 

J A 

dp t = d\” = —{ptAtPt + s £* s + *(Pt«t s - skjp t )}dt , 

where if T (s^,p^) = «(p£,s£). Using the integral form of the symmetric *-weakly 
continuous operators jl t , Vt '■ 0 —> S, one can obtain the stochastic integral repre¬ 
sentation of lnp t = f Q dAt(O) in Theorem 0 as well as the equation (15.711 and (15.811 
for tf t (Q = (£, At) and p t (£, £) = (£, p t £) with k t = p t + | s due to s T = -s, 

P tO ~ 2ikt{pt£, sO + i>M, s?) = ^(25?(£, (P‘ + \ s )Ct # . x )) 2 A(dx); 

25?(£, (Pt + \ s)C2, x ) = <6 2 P t3?Ct,x + s^Ct.x), Z e 5RS. 

Let’s point out that quantum filtering equations 15.71) . 15.81) . represented in the 
short form 


d$ t (£) + tiMtQdt = / 25ft(£, k t cS jX )dlt(dx) + u t (s£)dt, 

J A 

(5.9) -TTPt(£,£) + 2p t (£,a t £) = e t (s£,s£) + p, t (Pt£,Pt£), Po = 1, 

at 

where a* = ft t Pt + i&ts, give A*(£) = ln0 t (A £) in the form of the integral 

(5.10) X t (0 =i?o(| 0 ) + i fo + J^ {(s|, dT r ) + i (eM, s£) - M r (&■(£)> MO)) dr} 

over the stochastic trajectories = (p r (t,£) of the backward linear equation 15.41) 
with p r {£) =&r + Prt r , £o = £o(*> 0 and 

dT t =v t dt + S [ Ct x dYt(dx). 

J A 

Indeed, if d_^ r = £ r — £ r _ dr is the backward stochastic differential in (EH), then 
d{K,£ r ) = (Ir’M) + <i? r ,d_£ r ) and 

d(iv(£ r >tr)) = 2p r (C r ,d_| r ) +p r (| r ,t r .)dr. 

Using the equations EH and writting the equation 15.41) in the form d_^ r = 
o; ? .^ r d7- — 23?C,,dYV with respect to 

dF r (dx) = dU r (dx) + (p r £r, 23?Cr-,x)drA(dx) 
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- ^ - 2 

one can obtain by integrating by parts of the difference A*(£)— In p t — $o(£o) — 2 £o : 

^(0-^o(lo) + ^PtK,0-^S 

= {(l,di? + t p£dr) + (d_£,p(£))} 

= J {(4,dT) - (23?^!",^)+((a|,p|) + - jLt(p|,pO)dr} 

= L ^ s ^ ,dT ^ -^( 2K 0 d ^ + ^ (£(sis|)+/i(^,i?)-A(p(0 ; p(0))dr} 

= |( s C>dT) + t (e(s£,s£)-/S(p(£),p(t)))drj -lnp t , 

which gives 15.101) with lnp t = / 0 *(^(23fi£)dY — | p($, d)dr). 

Remark 4. Let us consider the case of the complex observation 12.1011 . Then 
the stochastic differential equation 

*d<M£) + {s(v u £) - (K t s£,d) + ^ e t (s£,s£)}0 t (£)}d t 
= 25? / <C t * lX ,^ t (0)dZ t (dx) + 29f / «(e,C t * lX )^(0d^(dx), 

JA J A 

corresponding to EB. has the Gaussian solution El) defined by the density 

p t = exp / {23? f ^ r (cf. x )d-^r-(dx) — £ T .(i? T .,'d r .)dr}, 

Jo J A 

where £*($,$) = f A |'i?(Ci. x ) | 2 ^(dx), and by the filtering equations 

d^t(C) + &t{iK t s£)dt = 23? / (£, x )dZ t (dx) + v t (s£)df, 

J A 

where dZ t (E) = dZ t (E) — dt L $ t (£ t x )A(dx). The corresponding complex form of 
Riccati equation is 
d 


df 


Pt(£,0 + 2p t (£, iKt s£) = £ t (s£, s£) - 2£ t (k t £, k^), 


where k t = p t + | s, kj = p t - | s. _ 

In the case of the complex Langevin equation pm . corresponding to iK t s = 
K t c on the invariant subspace S - , the posterior quasi-free dynamics with complex 
observation can be described in terms of the complex parameters n t = ^ £t. + i<^t , 
k\ = \ £ t - iU3 t , 

Xt = h |S", l t = kf|S"; Xt(0= 0, 1^ = 0, V£es+. 

In this case <j> t {£ © £ B ) = P t exp{i(x?(£) + Xt(£ # )) - ((*©*£>}, 

Pt = expj^ {23?^Xr(C!.,x) d ^r(dx) -£ r (x5.,Xr)}d?’ 

dXt(0 +Xt(«t<)di = /(!«?, Cb)d^(dx) + %«)dt, 

J A 

d t 1 , 1 x_ , 1 s 

— p t + p t « t c + c/«]p t = - ce t c- (p t - - c)e t (p t - - c), 
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where e t = Kt + k\ on 2 , p t = l t + \ c and 

dZ t (E) = dZ t (E) -At f x t (C t . x )A(dx). 

J E 

Example. Let us consider the stationary case £t = £, and a complete 

observation when A = B, £ t = £ is invertible and u> satisfies the c-conrmutativity 
condition luce = £cu> with e. One can easily prove that the last equation describes 
the continuous collapse of a posterior state to a Gaussian pure (coherent) state (j)^ 
with the minimal uncertainty Poo = y £ _1 |£c|, where 

£- 1 |ec| = e-V^cle- 1 / 2 = |c£|£-\ 6 = e 1 / 2 ^ 1 / 2 , |£| = (6 2 ) 1 / 2 . 

Indeed, if ujce = ecu;, then u>|c£| = |£c|uq because from the commutativity of c 
with £ _1 / 2 u;£ _ 1 / 2 there follows the commutativity of |c| with £~ 1 I 2 uj£~~ x I' 1 and 

W|C£| = o;£- 1 / 2 |£ 1 / 2 C £ 1 / 2 |£ 1 / 2 = E 1 / 2 ^ 1 / 2 ^ 1 / 2 ^- 1 / 2 ^ = |ec|«. 

Hence the Riccati equation 

d 1 

— Pt + l(ptWC - CUjpt) + Pf£Pt = - C£C, 

corresponding to e = £ — £ = 0, has the unique stationary positive solution p^ : 
Poo^c = c^Poo, Poo^Poo = \ £” 1/2 |c| 2 £~ 1 / 2 = \ c£c. The convergence p t -> p^ 
follows from properties of the Riccati equation with unique stationary solution 
Poo > 0. Thus in the case £ = el and uj = u> 1 the positive solution p t corresponding 
to po = 1 has the form 

|c| 1 + qt r , m 2 - |c| 

P t = — - , Qt = Qo exp{—e|c|t}, q 0 = , 

2 1 - qt 2 + |c| 

and Poo = \ |c|, Pt « Poo + |c|q 0 e _e|c|t for t > \ , if |c| > 0. Hence p t = Poo only 
in the purely quantum case |c| = 2, and p t —> 0 only in the purely classical case 
c = 0 when p t = 1/(1 + et). 

This result was obtained in jj for the case of positive definite c > 0 (a stable 
quantum system), when the stationary quantum linear filter 

dX t (0 + x t («x£)dt = [ (1£, Cx)( d ^t(dx) - Xt(CJMdx)dt), 

J A 

corresponding to ( t x = ( x , k = \ e + ilu, 1 = \ (|c| — c) and vt = 0, does not 
depend on the observable process Z f : 1 = 0, if c > 0. The a posterior state for the 
Gaussian initial wave function ijj 0 tends asymptotically to the ground state even 
without observation as the coherent a priori state: Xt(£) = Xo( e_KC *£) ~ 1 ' 0- 

In the contrary case c < 0, 1 = |c|, corresponding to the unstable system, the 
complete nondemolition observation is needed to keep the system in a state with 
the minimal uncertainty relation. This explains why the quantum open (unsta¬ 
ble) oscillator, corresponding j2J to the case S - = C = S + with |c| = 2, tends 
asymptotically to the pure Gaussian state under the continuous observation of its 
amplitude L = R (£), given by the measurement of the complex nondemolition 
process Z(t) = L(t) + W{t). 
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6. Conclusion 

The developed time-continuous quantum measurement theory based on the non¬ 
demolition principle as a superselection rule for the output observable processes 
abandons the von Neumann projection postulate as a redundancy in the stochastic 
extension of the quantum theory. It treats the wave packet reduction not as a real 
dynamical process but rather as the statistical inference of the posterior state de¬ 
scribed by the conditional wave-function for the prediction of the probabilities of 
the future measurements conditioned by the past observations. 

There is no need to postulate in quantum stochastic theory a nonunitary sto¬ 
chastic linear or nonlinear evolution for the continuous state-vector reduction as it 
is done in the phenomenological quantum theories of quantum trajectories, state 
diffusion or spontaneous localization 

puna mam nu. The nonunitary classical 
stochastic evolution giving the continuous reduction and localization of the poste¬ 
rior state can be rigorously derived mumm within the quantum stochastic unitary 
evolution of the correspondent compound system, the object of the measurement 
and the input Bose field in the vacuum state as it is shown here. 

References 

[1] Barchielli A. and Lupieri G. Quantum stochastic calculus, operation valued stochastic pro- 
cesses and continual measurement in quantum mechanics. Math. Phys., 26:2222-2230, 1985. 

[2] Belavkin V.P. Optimal measurements and control in quantum dynamical systems. Technical 
Report 411, Institute of Physics Copernicus University, Toruri, February 1979. 

[3] Belavkin V.P. Quantum filtering of Markovian signals with quantum white noises. Radiotech- 
nika i Electronica, 25:1445, 1980. 

[4] Belavkin V.P. Nondemolition measurements, nonlinear filtering and dynamic programming of 
quantum stochastic processes. In Modelling and Control of Systems, Proceedings Bellmann 
Continuous Workshop, Sophia-Antipolis 1988, Lect. Notes in Control and Inform. Sciences, 
vol.121,, pages 245—265. Blaquiere A., Springer-Verlag, Berlin-Heidelberg-New York-London- 
Paris-Tokyo, 1988. 

[5] Belavkin V.P. A new wave equation for a continuous nondemolition measurement. Phys. Lett. 
A, 140:355, 359, 1989. 

[6] Belavkin V.P. A posterior Schrodinger equation for continuous nondemolition measurement. 
Math. Phys., 1990. 

[7] Belavkin V.P. A quantum nonadapted Ito formula and generalized stochastic integration in 
Fock space. J. Funct. Anal., 1991. 

[8] Belavkin V.P. and A Barchielli. Measurements continuous in time and posteriori states in 
quantum mechanics. J. Phys. A. Math. Gen., 24:1495-1514, 1991. 

[9] Belavkin V.P. and Staszewski P. Watchdog effect for a free quantum particle. Phys. Rev. A., 
1991. 

[10] Diosi L. Continuous quantum measurement and Ito formalism. Phys. Lett. A., 129: 419-423, 
1988. 

[11] Ghirardi G.C., Pearle P., and Rimini A. Markov processes in Hilbert space and continuous 
spontaneous localization of systems of identical particles. Phys. Rev. A., 42: 78-89, 1990. 

[12] Bush, P., Lahti, P.J. Some remarks on unsharp quantum measurements, quantum nondemo¬ 
lition, and all that. Annal. Phys. 7, (5), 369-382 (1990) . 7, (5), 369-382 (1990) . 

[13] Bush, P.: Macroscopic quatum system and the objectification problem. In: Lahti, P., Mittel- 
staedt, P., (ed.). Symposium on Foundations of Modern Physics. Singapore: World Scientific 
1990. 

[14] Ghirardi G.C., A. Rimini, and T. Weber. Phys. Rev., 34D:470, 1986. 

[15] Hudson R.L. and Parthasarathy K.R. Quantum Ito’s formula and stochastic evolution, 93, 
301-323. Comm. Math. Phys., 1984. 

[16] Lindblad G. On the generators of quantum dynamical semigroups. Commun. math. Phys., 
48:119-130, 1976. 


QUANTUM CONTINUAL MEASUREMENT OF CCR 


25 


[17] Gizin N. Quantum measurements and stochastic processes. Phys. Rev. Lett., 52:1657-1660, 
1984. 

[18] Pearle P. Combining stochastic dynamical state-vector reduction with spontaneous localiza¬ 
tion. Phys. Rev. A., 39:2277-2289, 1989. 

Mathematics Department, University of Nottingham, UK 
E-mail address: vpbOmaths.nott .ac.uk 



